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SUPPRESSION OF CHEMOTACTIC EXPLOSION BY MIXING 


ALEXANDER KISELEV AND XIAOQIAN XU 


Abstract. Chemotaxis plays a crucial role in a variety of processes in biology and ecology. 
In many instances, processes involving chemical attraction take place in fluids. One of the 
most studied PDE models of chemotaxis is given by the Keller-Segel equation, which describes 
a population density of bacteria or mold which attract chemically to substance they secrete. 
Solutions of the Keller-Segel equation can exhibit dramatic collapsing behavior, where density 
concentrates positive mass in a measure zero region. A natural question is whether presence 
of fluid flow can affect singularity formation by mixing the bacteria thus making concentration 
harder to achieve. 

In this paper, we consider the parabolic-elliptic Keller-Segel equation in two and three di¬ 
mensions with additional advection term modeling ambient fluid flow. We prove that for any 
initial data, there exist incompressible fluid flows such that the solution to the equation stays 
globally regular. On the other hand, it is well known that when the fluid flow is absent, there 
exist initial data leading to finite time blow up. Thus presence of fluid flow can prevent the 
singularity formation. 

We discuss two classes of flows that have the explosion arresting property. Both classes 
are known as very efficient mixers. The first class are the relaxation enhancing (RE) flows 
of [5]. These flows are stationary. The second class of flows are the Yao-Zlatos near-optimal 
mixing flows |39j . which are time dependent. The proof is based on the nonlinear version of the 
relaxation enhancement construction of [5] , and on some variations of global regularity estimate 
for the Keller-Segel model. 


1. INTRODUCTION 

Chemotaxis is ubiquitous in biology and ecology. This term is used to describe motion 
where cells or species sense and attempt to move towards higher (or lower) concentration of 
some chemical. The hrst mathematically rigorous studies of chemotaxis effects have been by 
Patlak [3l] and Keller-Segel [19], [20]. The latter work involved derivation and hrst analysis 
of the Keller-Segel system, the most studied model of chemotaxis. The Keller-Segel equation 
describes a population of bacteria or mold that secrete a chemical and are attracted by it. In 
one version of the simplihed parabolic-elliptic form, this equation can be written in as (see 
e.g. m) 

(1.1) dtp-Ap + V ■ {pV{-A)-^p) = 0, p{x,0) = po{x). 

The last term on the left hand side describes attraction of p by the chemical whose concentration 
is given by c{x,t) = {—A)~^p{x,t). The literature on the Keller-Segel equation is enormous. It 
is known that in dimensions larger than one, solutions to (HU can concentrate hnite mass in 
a measure zero region and so blow up in hnite time. We refer to ra. n. and [12] for more 
details and further references. 

Typically, chemotactic processes take place in huid, and often the agents involved in chemo¬ 
taxis are also advected by the ambient how. Some of the examples involve monocytes using 
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chemokine signalling to concentrate at a source of infection (see e.g. IH113I]), sperm and eggs 
of marine animals practicing broadcast spawning in the ocean (see e.g. EES), and other nu¬ 
merous instances in biology and ecology. Our goal in this paper is to study the possible effects 
resulting from interaction of chemotactic and fluid transport processes. Of particular interest 
to us is the possibility of suppression of hnite time blow up due to the mixing effect of fluid 
flow. The problem of chemotaxis in fluid flow has been studied before; for example, in a setting 
similar to ours [21], [22] studied the effect of chemotaxis and fluid advection on the efficiency 
of absorbing reaction. Moreover, in a series of papers [2^, [TO] . i. IZSl. |2i a very interesting 
problem coupling chemotactic density with fluid mechanics equations actively forced by this 
density has been considered in a variety of different settings. The active coupling makes the 
system more challenging to analyze, but in some cases intriguing results involving global ex¬ 
istence of weak solutions (the dehnition of which implies lack of the h-function blow up) have 
been proved. These results, however, apply either in the setting where the initial data is small 
(see e.g. [28]) or close to constant ra. or in the systems where both chemotactic equation and 
the fluid equation have globally regular solutions if not coupled. In other words, to the best of 
our knowledge, there have been no rigorous results providing an example of suppression of the 
chemotactic explosion by fluid flow; only results showing that presence of fluid flow does not 
lead to blow up for the initial data that would not blow up without the flow. 

In this paper, our main focus will be on the question whether incompressible fluid flow can 
arrest the hnite time blow up phenomenon which is the key signature of the Keller-Segel model. 
There are two possible huid how ehects that can be helpful in hnite time blow up prevention. 
The hrst applies in inhnite regions, where strong how can help dihusion quickly spread the 
density so thin that chemotactic ehects become weak. The second ehect is more universal and 
subtle to analyze, and involves mixing in a hnite volume. In this case, the concentration may 
remain signihcant, but the how is constantly mixing density and preventing chemotaxis from 
building up a concentration peak. We are primarily interested in the mixing ehect, and so we 
will consider a hnite region setting. It will also be convenient for us to adopt periodic boundary 
conditions and to consider the Keller-Segel equation with advection on a torus. This is not 
essential, and many of our results also apply on a hnite region with Neumann, Dirichlet or 
Robin boundary conditions. 

Let us now state our main result. Since we are working on T'^, we will dehne the concentration 
of the chemical by factoring out a constant background: c{x,t) = {—A)~^{p{x,t) — p). Here 
p{x, t) G is the species density, and p is its mean over T'^. The inverse Laplacian can 
be dehned on the Fourier side, or by an appropriate convolution as will be discussed below. 
Consider the equation 

(1.2) dtp + {u-V)p - Ap+ V ■ {pV{-A)~^{p - p)) = Q, p{x,Q) = pt){x), x e 
We will prove the following theorem. 

Theorem 1.1. Given any initial data po 0, po E C°°(T‘^), d = 2 or 3, there exist smooth 
incompressible flows u such that the unique solution p{x, t) of / If. ill) is globally regular in time. 

Remark 1.1. 1. We will provide more details on the specihc choice of the hows later. 

2. The restriction d = 2 or 3 is essential for our method to work. The case d = 4 is in some 
sense borderline for our estimates: the needed estimates just fail. Extending our approach to 
this case, and especially to d > 4, would require new ideas. 
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Please check Remark 12.11 after Theorem 12.11 for more discussion. 

3. As we already mentioned, the periodic boundary condition is not very essential. One can 
get a similar result in the bounded domain; one issue is that examples of stationary RE flows 
are not available in a general bounded domain. On the other hand, the Yao-Zlatos flows can 
be constructed in general bounded domains and so can be used to achieve an analogous result 
in this more general setting. See Remark 16.31 for more details. 

4. The theorem certainly holds under weaker assumptions on the regularity of the initial data. 
In this paper, for the sake of presentation, we do not make an effort to optimize the regularity 
conditions. The scheme of our proof and the kinds of the flow examples that we have will make 
the connection between mixing properties of the flow and its ability to suppress the chemotactic 
blow up quite explicit. 

It is well known that the solutions to the Keller-Segel equation can form singularities in hnite 
time. The hrst rigorous proof of this result in the case where domain is a two-dimensional disk 
was given by Jager and Luckhaus [IB]. Their proof is based on radial geometry and comparison 
principles. Nagai [33] has provided a proof of hnite time blow up in more general bounded 
domains. We have not found a hnite time blow up proof for the periodic case in the literature. 
Although it can be obtained by modihcation of the existing arguments, in the appendix we will 
provide a short independent construction of such examples in the case of two spatial dimensions. 
This will imply that Theorem 11.11 indeed provides examples of the suppression of chemotactic 
explosion by huid mixing. 

We note that huid advection has been conjectured to regularize singular nonlinear dynamics 
before. The most notable example is the case of the 3D Navier-Stokes and Euler equations. 
Constantin [1] has proved possibility of hnite time blow up for the 3D Euler equation in if 
the pure advection term in the vorticity formulation is removed from the equation. Hon and 
Lei have obtained numerical evidence for the hnite time blow up in a system obtained from the 
3D Navier-Stokes equation by the removal of the pure transport terms [13]. In fact, hnite time 
blow up has been also proved rigorously in some related modihed model settings [15], [16], [14] . 
Of course, the proof of the global regularity of the 3D Navier-Stokes remains an outstanding 
open problem, so whether the 3D Navier-Stokes equation exhibits “advection regularization” is 
an open question. See also [21] for more discussion. As another example of related philosophy, 
we mention the paper [2] on the elliptic problem with “explosion” type reaction. There is 
no time variable and so no hnite time blow up in this paper, but it shows that certain hows 
can signihcantly ahect the “explosion threshold”: namely, the value of the reaction coupling 
parameter beyond which there exist no regular positive solutions. To the best of our knowledge, 
the examples that we construct here are the hrst rigorous examples of the suppression of hnite 
time blow up by huid mixing in nonlinear evolution setting. It should be possible to extend 
our method to cover some other situations, which we will briehy discuss in the last section. 

The paper is organized as follows. In Section 2 we prove the -based global regularity 
criterion that we will use. In particular, it guarantees global regularity as far as the norm 
of the solutions remains bounded. In Section 3, we set up the strategy for controlling norm 
via norm of the solution. Despite the presence of nonlinearity, due to the dissipation term, 
if the norm is sufhciently large, the norm has to decay. Basically, what this implies is 
that the solution to the Keller-Segel equation cannot blow up in the manner which makes 
norm of the solution grow much faster than the norm. In Section 4, we prove the key result 
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on approximation of the solution of the Keller-Segel equation by solution of pure advection 
equation on small time scales. This result sets up our strategy to use the fluid advection to 
drive up the norm of the solution to the full Keller-Segel equation by fluid mixing in order to 
keep the norm in check. Indeed, the highly mixing flows are known to increase the norm 
of the solution of the advected passive scalar. If this solution stays close to the solution of the 
nonlinear problem long enough, this guarantees that the norm of the nonlinear solution will 
be large, too. In Section 5, we put together all the components of the argument and prove that 
the relaxation enhancing flows of [5] suppress chemotactic blow up. We focus on the case of 
weakly mixing flows. In Section 6 we outline another example of flows with analogous property, 
the Yao-Zlatos efficient mixing flows. Finally, in Section 7 we briefly discuss possible future 
extensions. There are two appendices: the hrst one contains construction of hnite time blow 
up in the Keller-Segel equation with periodic data, and the second includes sketches of proofs 
of some functional inequalities that we use in the paper. 

Throughout the paper, C will stand for universal constants that may change from line to 
line. 

Acknowledgment. AK has been partially supported by the NSF grant DMS-1412023. 
XX has been partially supported by the NSF grant DMS-1535653. We thank Eitan Tadmor, 
Changhui Tan, Yao Yao and Andrej Zlatos for helpful discussions. 

2. Global existence: the criterion 

In this section, we will show that to get the global regularity of fll.2p . we only need to have 
certain control of the spatial norm of the solution. The following theorem is a direct analog 
of Theorem 3.1 in [22], where it was proved in the setting. We will provide a sketch of the 
proof for the sake of completeness. Throughout the paper, we will use notation for the 
homogeneous Sobolev space in spatial coordinates, that is we set 

ii/iih= X 

Theorem 2.1. Suppose that pQ G d = 2 or 3, po ^ 0, and suppose that u G C°° is 

divergence free. Assume [0,T] is the maximal interval of existence of the unique smooth solution 
p{x,t) of the equation M.2t) . Then we must have 

(2.1) ^ ||p(-, r) - pII oo. 

Remark 2.1. In other words, the smooth solution can be continued as far as the integral in 
time of the appropriate power of the Lf norm of solution in space stays hnite. Note that the 
mean value of p is conserved by the evolution, so p(-,t) = Po- We will denote it p throughout 
the rest of the paper. One may or may not include p into fl2.ll) . these criteria are equivalent. 
One can verify that the scaling of fl2.ip is sharp in the sense that it is a critical quantity for 
fll.2l) . One way to see this criticality is through an informal scaling argument or dimensional 
analysis. Indeed, one can check that A^p(|, is a solution to fll.2p (with a different period) for 
every A as long as p(a;, t) is a solution. The quantity in fl2.ip with respect to this scaling. This 
observation also means that for d ^ 4 the boundedness of the norm may not be sufficient 
in order to get the regularity of the solution to fll.21) . 
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Proof. The existence and nniqneness of smooth local solntion can be proved by standard meth¬ 
ods, so we will focns on global regnlarity. Let s ^ 2 be integer. Mnltiply fll.2p by (—A)^p and 
integrate. We get 


-dt 

2 




< 


Ijd 


(Vp).(V(-A)-'(p-p))(-A)><ij^ 


Ijd 


pip- p)i-^ypdx 


+ C\\u\ 




2 

ffs+l- 


Here we integrated by parts s times and used incompressibility of u to obtain 



{u ■ V)pi—Aypdx 


< 


C'l|w||cH|p|| 


2 

H=' 


Consider the expression 



pip - p)i-^ypdx. 


Integrating by parts, we can represent this integral as a sum of terms of the form 


where I 



D^pD^-\p- p)D^pdx, 


0,1,..., s and D denotes any partial derivative. By Holder inequality, we have 


[ D‘pD‘-\p - p)D-pdx \\D‘p\\„\\D‘-‘(p - p)\\„\\p\\f,., 

with any 2 ^ pi, qi ^ oo satisfying = 1/2. For any integer 0 ^ m ^ n, and mean zero 

/ G C'°°(T^), we have Gagliardo-Nirenberg inequality 


( 2 , 2 ) 




a 

h"’ 


m 


a = 


d I d 

p ^ 2 


n 


which holds for 2 ^ p ^ oo unless a = 1, and if a = 1 for 2 ^ p < oo. We will sketch a short 
proof of fl2.2l) in the appendix to make the paper more self-contained. 

Take = y, O’/ = ^- Then for all I > 0, applying (12.2p . we get 


2s —<i+4 2s+d 

\\D^p\M\D^-yp-mLn^c\\p-p\\it^^ iipiig::f. 

In the I = 0 case, we use 

IIpIU- ^ Up - pIU- + p ^ c||p - p\\l^ l|p|||:+l^ + p- 

Therefore 

pip-p)i-^ypdx 

Next, consider 

[ (Vp)-(V(-A)-i(p-p))(-Arpdx. 



IIp-pII 


1 - 

L2 


2(s+l) 


d 

2(s+l) 

ps + 1 


+ P 


S +1 

p-s+1 


iip-pii/r 


Hs- 
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Integrating by parts s times, we get terms that can be estimated similarly to the previous case, 
using the fact that the double Riesz transform dij{—A)~^ is bounded on L^, 1 < p < oo. The 
only exceptional terms that appear which have different structure are 


Jrd 

But these can be reduced to 


{di,...diyp) ■ (V(-A) {p- p))di^...di^pdx. 


Ijd 


{di^...di,pf{p- p)dx 


by another integration by parts, and estimated as before. Altogether, we get 


(2.3) 


-dt 

2 




1 - 


+ C\\u\\c4p\\%s 


2(s+l) 


2(s+l) 

Lfs+i 


+ P 


S+l 


\\P-P\\ir 


2 

ijs+l- 




Observe that 

(2.4) 

Split the hrst term on the right hand side of (12.31) into two parts, and estimate them as follows. 
First, 

c||p-p|l2^l|p||jSflLIL. ^ l|p-/5|U.||p||y„|Ll47 

«jiipiiy«+ci[p-piii?iipiiy. 

Second, 

p\\p\\ps+i\\p- pWIVWpWh^ ^ Pllp|lir»l|p|lir=+i 

siiipiiy«+c7iipfA.- 

We used the Poincare inequality and (12.41) in the hrst step. Recall the following Nash-type 
inequality 

(2.5) IIpIIr. ^qipll^llplll^, 

the proof of which will be sketched in the appendix. Since p{x,t) ^ 0 and hence ||p(-,t)||Li = 
p > 0 is conserved in time, putting all estimates into fl2.3p we get 


( 2 . 6 ) 


:d, 


^ C 


Up - p|ll2‘' + p^ + 


\u\\c‘ 


- cp 


2s + d 


2s+d 


From this differential inequality and integrability of ||p(-,t) — p\\'^ 2 'ljd^ in time, a hnite upper 
bound for ||p(-,t)||Lf 3 follows for all times. In fact, due to the last term on the right hand side 
of fl2.6p . it is not hard to show that there is a global, not growing in time, upper bound for any 
norm of p. □ 
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3. An condition for an absorbing set in 


Due to Theorem 12.11 to show global regularity of the solution p{x,t) to (IIAh . it suffices to 
control its norm in spatial variables. In this section, we prove a simple criterion that says 
that if the norm of a solution is sufficiently large compared to its norm, then in fact the 

norm is decaying. Our overall strategy will be then to show that mixing can increase and 
sustain the norm of solution. This will block the norm from ever growing too much, 
leading to global regularity. 

Proposition 3.1. Let p{x,t) ^ 0 be smooth local solution to U.W set on T'^, d = 2 or 3. 
Suppose that ||p(-,t) — p||l 2 = B > 0 for some t^O. Then there exists a universal constant Ci 
such that 

(3.1) ||p(-, t + t) — p\\i ,2 ^ 2B for every 0 ^ r ^ Oi min(l, B~^^). 

Moreover, there exists a universal constant Cq such that if, in addition, 

(3.2) \\p{;t)\\%, ^ Bl ^ CoB^-^ + 2pB\ 
then dt\\p{-,t)\\L 2 < 0. 

Remark 3.1. 1. In particular, due to Theorem 12.11 it follows that if ||p(-,t) — p\\l'^ ^ B, then 
the local smooth solution persists at least till t + Ci min ^1, p~^, 

2. Here and below, by universal constant we mean a constant that does not depend on any 
parameters of the problem: namely, a constant independent of u, Pq, and B. 


Proof. Let us multiply both sides of (II.2p hy p — p and integrate. Then 
1 
2 

Observe that 


(3.3) ndt\\p- p\\l 2 =-\\p\\^r+ I V-(pV(-A) {p-p)){p-p)dx. 


ijd 


ijd 


pV(—A) {p — p)Wpdx= / p{p — p)dx— Vp-V(—A) {p — p)dx 


I Yd 


I Yd 


Therefore, 


I Yd 


pV (-A) {p- p)Vpdx = - p {p-p)dx, 


I Yd 


and then the integral on the right hand side of fl3.3p is eqnal to 


- V-(pV(-A) {p-p){p-p))dx=- p{p-p)dx, 


I Yd 


I Yd 


Next, notice that 


p^{p — p)dx= I {p — p)^dx + 2p I {p — pYdx — 2p^. 


I Yd 


I Yd 


I Yd 


By a Gagliardo-Nirenberg ineqnality (see e.g. [3T] or [22] for a simple proof), we have 

n d d 12 — 2 id 

llp-p|li3 ^ c\\p-p\\~ 2 ^\\p\\i, ^ iipII^i + c'iiip-p||/T'^, 
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where in the second step we applied Young’s inequality. Applying all these estimates to fl3.3p 
yields 

12-2d 

(3-4) p||^2 ^-||p||^i + CqIIp - pll/a"'* + 2p||p - p||i2. 

Solving the differential equation 

/'(r) = C’/^(r) + 2p/(r), /(0)=i?^ 


leads to the solution 
(3.5) 


/M 


B'^ exp(2pr) 


A standard comparison argument can be used to show that ||p(-, t + t) — p |||2 ^ /('t)- On the 
other hand, a straightforward estimate using fl3.5p gives existence of a constant Ci such that if 
r ^ Cl min(l, p~^, then /(r) ^ iB^. 

The second statement of the lemma follows directly from fl3.4p and an assumption po ^ 0 
(which implies p ^ 0). 

□ 


4. An approximation lemma 

We can now outline our general strategy of the proof of chemotactic blow up suppression in 
more detail. We know that the control of norm in spatial coordinates is sufficient for global 
regularity. We also see that if norm of the solution is large, then norm is not growing. 
On the other hand, flows with strong mixing properties tend to increase norm of solution. 
Hence our plan will be to deploy such flows, at a sufficiently strong intensity, to make sure 
that the norm of the solution stays high, at least whenever the norm is not small. The 
hrst hurdle we face, however, is to show that the mixing property of flow persists in the full 
nonlinear Keller-Segel equation. 

In this section we prove a key result on the approximation of solutions to the Keller-Segel 
equation with advection fll.2p by solutions of the pure advection equation. We will be looking at 
the intense advection regime, and consider small, relative to all parameters except the strength 
of advection, time intervals. It is natural to assume that in this case most of the dynamics we 
observe is due to advection, though the exact statement of the result requires care since both 
diffusion and chemotactic terms are not trivial perturbations. 

Let us consider the equation fll.21) 

dtp + (u ■ V)p - Ap + V ■ (pV(-A)“^(p - p)) = 0, p(x, 0) = po(x), 

X G with d = 2,3. We will assume that the vector held u is divergence free and Lipschitz 
in spatial variables. It may be stationary or time dependent. Let us denote p(x, t) the unique 
Lipschitz solution of the equation 

(4.1) dtri + {u-V)r] = 0, p(x, 0) = po(x). 

If we dehne the trajectories map by 

(4.2) ^ $o(a;) = x, 
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then r]{x,t) = po($t 

We start with the following simple observation. Denote the Lipschitz semi-norm 


Wfhip = sup 

^,y 


k - y\ 


Lemma 4.1. Suppose that the vector field u is incompressible and Lipschitz in spatial variable 
for each t ^ 0, \\u{-,t)\\np(jd^ ^ D(t), D(t) G oo). Let ri{x,t) be the solution of ( (^. j[ j. 

Then for every t ^ 0, and for every pq G H^, we have 

(4.3) llh(L^)lliti ^ ^(^)IIPoIIhi, where F{t) = exp D{s)ds 

Proof. If u is incompressible and Lipschitz in spatial variable for each time, then the trajectories 
map is area preserving, Lipschitz in x and invertible for each t, and the inverse map 
is also Lipschitz in spatial variables. Moreover, ||*hL^IUjp ^ exp(C' Jq D{s)ds) (see e.g. [3U]L 
The evolution p{x,t) = is a Lipschitz regular coordinate change of an function 

Pq. The bound fl4.3l) follows from the well known properties of functions under Lipschitz 
transformations of coordinates [iQ]. □ 



We are now ready to prove the approximation lemma. 

Lemma 4.2. Suppose that the vector field u{x,t) is incompressible and Lipschitz in x, and is 
such that (^7^ is satisfied with F{t) G L“^[0,oo). Let p{x,t), p{x,t) be the solutions of / li.gj) . 
63 respectively with po ^ 0 G Lfh Suppose that the unique local smooth solution p{x, t) exists 
for t G [0, T]. Then for every t G [0, T] we have 

(4.4) L\\p - nWh < -IIpII?,. + 4||p„||? + Clip - Plli, (llp - Pll* + p^j . 

Of course a direct analog of lemma holds without assumption po ^ 0; we only need to replace 
p2 in (ITTD with IIpII^i. 

Proof. A direct computation using divergence free property of u shows that 

= / Ap{p-p)dx- [ V ■ {pV{-A)-\p-p)){p-p)dx 

< -IIpIIL + llpIlH.lkllp. + l|pV(-A)-‘(p - p)||„||p||p, 

+ \\p^{-^)~\p-P)h4v\\m- 

Applying Holder and Gagliardo-Nirenberg inequalities, we can estimate 
||pV(-A)-^(p-p)|U= ^ ||p|U3||V(-A)-^(p-p)|Ue 

^ (IIpIiIjIp - pW'l'J + p) IIV(-A)-^(p- p)|||j|V(-A)-i(p- p)C 
< c (l|p|||il|p-p|ll2® +p) ||p-p||l2- 

Here the last step follows from simple estimates on Fourier side. Given these estimates, several 
applications of Young’s inequality show that the right hand side of fl4.5p can be bounded above 


"SIM 
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by 

- ibiiL+ jIIpiil + iipiiy+giipiiy+ c\\p - 
+ gllpliy + c'lip — p|ii2p^ + iipUL- 

With help of Lemma 14.11 the estimate fl4.4p quickly follows. □ 

5. Proof of the main theorem: the relaxation enhancing flows 

Our first example of flows that can stop chemotactic explosion will be relaxation enhancing 
flows of [5] . These stationary in time flows have been shown to be very efficient in speeding up 
convergence to the mean of the solutions of diffusion-advection equation. A particular type of 
the RE flows are weakly mixing flows, a well known class in dynamical systems theory which 
is intermediate in mixing properties between mixing and ergodic jS]. Let us briefly review the 
relevant definitions. 

Given an incompressible vector field u{x) which is Lipschitz in spatial variables, recall defini¬ 
tion fl4.2p for the trajectories map Then define a unitary operator U^f{x) = /(<h^^(a:)) 

acting on L^(T'^). 

Definition 5.1. The flow u{x) is called weakly mixing if the spectrum of the operator U = 
is purely continuous. 

The flow u{x) is called relaxation enhancing (RE) if the operator U (or properly defined 
operator (n • V)) has no eigenfunctions in other than a constant function. 

Remark. The fact that we talk about the spectrum of U rather than (u ■ V) is a minor technical 
point. The symmetric operator i{u-V) is unbounded on L^, and sometimes needs to be extended 
to its natural domain (rather than just H^) to become self-adjoint and to be a true generator 
for U. To avoid these technicalities, it is convenient to make the definition in terms of U, which 
is bounded and so can be defined on smooth functions and then extended to the entire by 
continuity. 

Examples of weakly mixing flows on are classical and go back to von Neumann |3B] (just 
continuous u{x)) and Kolmogorov [23] (smooth u{x)). The Kolmogorov construction is based 
on the irrational rotation on the torus with appropriately selected invariant measure. Lack 
of eigenfunctions is established by analysis of a small denominator problem on the Fourier 
side bearing some similarity to the core of the KAM theory. The original examples are not 
incompressible with respect to the Lebesgue measure on but a smooth change of coordinates 
can be applied to obtain incompressible flows with the same properties. Weakly mixing flows 
are also RE, but there are smooth RE flows which are not weakly mixing: these do have 
eigenfunctions but rough ones, lying in but not in (see [5] for more details and examples). 

We now state our first main theorem. Consider the equation 

(5.1) A(n • V)p^ - -F V ■ (p"^V(-A)"^(p - p)) = 0, p"^(a:, 0) = po(a;). 

Here A is the coupling constant regulating strength of the fluid flow that we will assume to be 
large. We note that dividing the equation by A and changing time, we can instead think of 
all the results below as applicable in the regime of weak diffusion and chemotaxis on long time 
scales. 
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Theorem 5.2. Suppose that u is smooth and incompressible vector field on d = 2,3, which 
is also relaxation enhancing. Assume that p ^ 0 G C'°°(T“'). Then there exists an amplitude Aq 
which depends only on po ond u such that for every A^ Aq the solution p^[x,t) of the equation 
L5.1\) is globally regular. 

We will only prove Theorem 15.21 in the case of weakly mixing flows. This serves our main 
purpose of providing an example of chemotactic blow up-arresting flow. In the general RE case, 
the proof is a fairly straightforward extension of an argument for the weakly mixing flow and 
the point spectrum estimates in [Sj (namely, Lemma 3.3 and part of the proof of Theorem 1.4 
dealing with point spectrum). 

Before starting the proof, we need one auxiliary result from [5]. Let Pn be the orthogonal 
projection operator on the subspace formed by Fourier modes \k\ ^ N: 

PNf{x) = 

\k\^N 

Lemma 5.3. Let U be a unitary operator with purely continuous spectrum defined on L^(T“'). 
Let S = {(j) E L"^ : ||0||l2 = 1}, and let K <Z S be a compact set. Then for every N and every 
a > 0, there exists Tc{N, a, K, U) such that for all T ^ Tc{N, a, K, U) and every f E K, we 
have 

(5.2) WPNU^ffdt^a. 

This lemma connects the issues we are studying with one of the themes in quantum mechanics, 
namely the propagation rate of wave packets corresponding to continuous spectrum. Lemma 
15.31 is an extension of the well-known RAGE theorem (see e.g. H) which is a rigorous variant of 
a folklore quantum mechanics statement that quantum states corresponding to the continuous 
spectrum travel to infinity. In our case, travel to infinity happens not in physical space, but 
in the modes of the operator —A (that is, in Fourier modes). We refer to j3] for the proof of 
Lemma 15.31 

Now we are ready to give the proof of Theorem 15.21 

Proof of Theorem \5.^ Fix any B > ||po — Po||l 2 . If for all times we have that \\p^{-,t) — p\\l^ < 
B, then the solution stays globally regular by Theorem 12.11 Otherwise, let 

to = inf{t| ||p^(-, t) - p\\l 2 = B}. 

Since the solution is smooth up to time to, we also have that ||p^(-,to) — p||l 2 = B; thus to is 
the first time the norm of p^{x,t) — p reaches B. Note that by Proposition 13.11 we must 
also have ||p"^(-, to)||_ffi < Bi, where Bi = CqB -|- 2pB‘^. 

We are going to show that if A ^ Aq{B,p,u) is sufficiently large, then after a small time 
interval of length r that we will define shortly, we will have ||p"^(-, to + u) — p||l 2 < B. Moreover, 
we will have ||p"^(-, t) — p \\]^2 ^ 2B for every t E [to, to + r]. This will prove the theorem, as the 
argument can be applied repeatedly each time the L^ norm reaches the level B, showing that 
Wp'^i-A) - pWl^ ^2B for all times. 

Denote the eigenvalues of —A on in an increasing order, 0 = Ai ^ A 2 ^ ... ^ A„ ^ ... 
Choose N so that 

(5.3) \n > 16Go(2R)ra + 32p, 













12 


ALEXANDER KISELEV AND XIAOQIAN XU 


where Cq is the constant appearing in fl3.4p . Observe that XnB^ > Bf. Dehne the compact set 
it: C ^ by 

K = {(j) e S\\m%, ^ X^} 

(recall S is the unit sphere in L^). Let U be the unitary operator associated with our weakly 
mixing flow u as above. Fix a = 0.01. Let Tc{N,a, K,U) be the time threshold provided by 
Lemma 15.31 

We proceed to impose the first condition on Ao{po, u). We define r as below and require that 


(5.4) 


^ ( AT 


for every A ^ Ao, where Ci is the constant appearing in Proposition 13.11 in fl3.ip . It follows 
from Proposition 13.11 and Theorem 12.II that \\p^l-,t) — p)\\l 2 < 25 for t G [to,to + t] and so p^ 
remains smooth on this time interval. 

Let us introduce a short-cut notation (po^x) = p^{x,to). Let p^{x,t) be the solution of the 
equation 

dtp^ + A{u-V)p^ = 0, ?7^(x, O) = 0o- 

Then p^{x,t) = U^^(f)o, and we have 


(5.5) -[ ||PAr?7^(x,t)||^2dt = - / \\PNU^^(po\\l2dt = ^ [ \\PNU"(po\\hds ^ aB‘^. 

Jo Jo Jo 

Here we applied Lemma 1531 to the vector 0 o/||0o||l 2- Indeed, 0 o/||0o||l2 G K since by fl5.3p we 
have 

|| 0 o|||-i ^ Pi < < AAr|l0o||i2. 

Note also that fl5.4l) ensures the applicability of Lemma 15.31 and the validity of the last bound 
in (15.5p . 

We now impose the second and last condition on Aq. It will be convenient for us now to 
denote by t time elapsed since to- By fhe approximation Lemma [4.21 and since we know that 
||p'^(-,to)|||.i ^ Pi < AatP^, as well as ||p^(-, to-h t) - p\\l^ ^ 25, we have 


(5.6) ^0 + t)- v%, t)\\l, ^ AXNB‘^F{Atf + ^5^(5^ + p^) 
for all t G [0,r]. Here r = Tc/A as before. Choose Ho so that 

(5.7) ^ ° F{s)^ds + CriB^d + p^) < 0.01 

A Jo 

for every A ^ Ag. Note that since u is smooth, F(t) is a locally bounded function. 

We claim that if H ^ Ho, then ||p^(-, to -l- r) — p||x ,2 ^5. First, the condition (15.7p allows us 
to control ||p"^(-,to -l- t) — p||r 2 more tightly, which is convenient. Indeed, since ||? 7 '^(-, t)|| r 2 = 
|10o||l 2 = 5 for all t > 0, ([521) and (EIP imply that 


(5.8) 


0.95 ^ ||p"^(-, to -l- t)||i2 ^ 1.15 


for t G [0,r]. Furthermore, by the estimates (15.5p . (15.Op and (15.7p we have 

-/ \\PNp^{-,to + t)\\l2dt [ WPNv'^i-A^hdt 
'T Jo 'T Jo 
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Combining this estimate with fIS.Sp we obtain 

(5.9) - f ||p"^(-,to + ^ - [ ^N\\{d — PN)p^{-,to + t)\\‘L'2dt ^ -XnB'^- 

Jo 'T Jo ^ 

Now we come back to (13.41) : 

12-2d 

dt\\p^-p\\l2 ^ -IIp^II^i +<^011^^ -p||/2"‘' +2p\\p^ -P\\l2. 

By the estimates (15.91) . (15.Sp and (15.3p . we have 


(5.10) 


^0 + u) ~ p|li2 ^ + T + Cq{2B) 4-d + 2p{2B)^'^ 


< 



B^ ^ B^. 


This completes the proof. 


□ 


We see that the bound we obtained on the decay of the norm in (I5.10p is stronger than what 
we needed. In fact, with slightly more effort we can obtain stronger results. We now present 
an extension of Theorem 15.21 that establishes a complete analog of “relaxation enhancement” 
established in [5] for the diffusion-advection equation for the case that also includes chemotaxis. 
Namely, we show that not only fluid flow can prevent hnite time blow up, but in fact it can 
enforce convergence of the solution to its mean in the long time limit. Intense fluid flow can 
also act to create an arbitrary strong and fast drop of \\p^{-,t) — P\\l 2 . 

Theorem 5.4. Suppose 0 ^ po E C°°(T'^), and let p^{x,t) be the solution of the equation H5.1\) . 
Let u be smooth, incompressible, relaxation enhancing flow. If Aq{po,u) is the threshold value 
defined in the proof of Theorem \5.A then for every A ^ Aq, we have 

(5.11) \\p^{-,t) - P\\l2 ^ 0 

as t —)■ cxo. The convergence rate is exponential in time, and can be made arbitrary fast by 
increasing the value of A. Namely, for every h > 0 and k> 0, there exists Ai = Ai{pq,u, k, 6) 
such that if A ^ Ai, then 

(5.12) ||p^(-, t) - p\\l2 ^ IIpo - 
for all t ^ 6. 

Proof. Let us prove (I5.12p . since (15.111) follows from similar (and easier) arguments. The proof 
largely follows the argument in the proof of Theorem 15.21 but let us outline the necessary 
adjustments. We set Bq = ||po — p\\l^- Choose N so that 

(5.13) Aat ^ max (^IOOk, 16Co(25o)^ + 32p) . 

Dehne the set K, as before, by {0 G S'|||0||^^ ^ Atv}. 
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For all times t where p^{x,t)/\\p^{-,t)\\L 2 ^ K, we have ||p^(-, t)||^i ^ XN\\p'^{-,t) - p\\l^. It 
follows from fl3.4p and 05.131) that at such times we have 

12-2d 

9t\\p - p\\‘l2 ^+ CoWp"^ - p\\j\~‘^ +2p\\p^ - p\\l2 
(5.14) . 

^ Av + <^0(2^0)+ 2 pj IIp"^ — p ||^2 ^ — pllia- 

Here in the second inequality we used that \\p^{-,t) — p\\i 2 ^ 2Bo for all times, as we know 
from the proof of Theorem 15.21 we can ensure by making A sufficiently large; we note that 
this bound will also follow from our argument below. Thus on the time intervals where 
p^{x,t)/\\p^{x,t)\\L^ ^ K we have exponential decay of ||p"^(-,f) — p\\l 2 at rate that would 
imply 05.12P if all times were like that. 

Suppose now that to is the smallest time such that p^{x,to) ^ K {to could equal 0). Let 
Tc{N,a, K,U) be the time threshold provided by Lemma 15731 (we set a = 0.01 as before). 
Repeat all the steps in the proof of Theorem 15.21 from dehning the time step r fl5.4p to 05.101) . 
with B replaced by ||p"^(-,to) — pWl'^- In addition, require that A is large enough so that 


(5.15) 

We arrive at the estimate 


T,{N,a,K, U) 
A 


^ (5/2. 


(5.16) 


||p"^(-,^o + t) - PIU2 ^ 


(^1 - ||p"^(-,io) -pII 


2 

L2 




Note that even though we do not control the norm of the solution for all times t, from 03.4p 
and 05.13P it is clear that for every t E [to,to + x] = Iq, we have 

(5-17) IIp'‘(V)-pHL « 

We continue further in time in a similar fashion. If p^{x,t)/\\p^{-,t)\\L^ ^ K, we have 05.141) . 
On the other hand, if 

tn = inf{f|f > tn-i + t,p^{xA)/\\p^{-A)\\l^ e JF}, 

we can apply Lemma 15731 and Lemma W?2\ on + r] obtaining 

\\p^{-.tn + r)- p\\l2 ^ e--i^^^\\p^{-An) - p\\l2, 

(5.18) , 

IIp^(l^)-pIIl2 ^ - pII^s, for every te [tnAn + r], 

Now given any t ^ 6, we can represent 


[o,i] = »'u(ur.„/,). 


where W is the set of times in [0,f] outside all Note that 05.14p holds for every s G W. 
Combining 05.14p and 05.181) . we infer that for every t ^ 6, we have 


l|p''(V)-p|lL<e-I^»'‘A)e4^A||p„_p|| 


\h 


^ e 


— ^ Awt 


IIpo - 


pII 


2 

L2- 


This proves 05.12p . 


□ 
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6. Yao-Zlatos flow 


In this section, we describe another class of flows that are capable of suppressing the chemo- 
tactic explosion. These flows arise as examples of almost perfect mixers satisfying some sort 
of natural constraints. They have the advantage of being somewhat more explicitly dehned 
than the RE flows which may be harder to picture. For this reason we will be able to get an 
explicit estimate, albeit rather weak, for the intensity of mixing necessary to arrest the blow 
up as a function of the norm of the initial data. For the RE flows, such estimate would be 
difficult to obtain, primarily due to the challenge of estimating the time Tc from Lemma 15.31 
A quantitative estimate on Tc would require delicate spectral analysis of the operator u • V, 
something that for the moment is out of reach. On the other hand, in contrast to the RE flows, 
Yao-Zlatos flows are time dependent and active - their construction depends on the density 
being mixed. We remark that a re lated class of efficient mixer ffows has been also considered 
in [I]. 

We refer to [39] for a detailed discussion of different notions of mixing and the general 
background, and for further references. For our purpose here, we need one particular result 
from [39], that we set about to explain. Let = [—1/2,1/2)^. Consider the dyadic partition 
of with squares Qnij given by 


( 6 . 1 ) 


Qnij 


i + 1 


2 "’ 2 ^ 


X 


J_ j + 1 

2n’ 2" 


hj = -2 


n—1 


2 "-^ - 1 . 


Suppose /o G and is mean zero, and u is an incompressible flow which is Lipschitz 

regular in spatial variables. Let f{x,t) denote the solution of transport equation 


( 6 . 2 ) 


dtf + {u-V)f = 0, f{x,0) = fo{x). 


Theorem 6.1. [Yao-Zlatos] Given any a, e G (0,1/2], there exists an ineompressihle flow u 
such that the following holds. 

(6.3) 1. IIVm(-, f)||Loo ^ 1, for every t. 

2. Let n = [I log 2 (Ke)|] -|- 2, where [x] denotes the integer part of x. Then for some 

and every Qnij as in h6.1\) we have 


(6.4) 


\Q 


nij\ JQr, 


f{x,Tn,e) dx 


< «ll/( 


0 L° 


This theorem provides a flow u that satishes uniform in time Lipschitz constraint fl6.3p and 
mixes the initial density /o to scale e with the error at most k in time The construction 
in [59] employs a multi-scale cellular flow. On the nth stage of the construction, the goal is 
to make the mean of the function on each of Qnij close to zero, starting with n = 1. This 
is achieved by cellular flows which are designed to have the same rotation time period on 
streamlines away from a thin boundary layer. Such an arrangement makes evolution of density 
in each cell amenable to fairly precise control, and makes it possible to ensure that the “nearly 
mean zero” property of the solution propagates to smaller and smaller scales. We refer to [39] 
for the details. Below, we outline only some adjustments that are needed to obtain Theorem 16.II 
from the arguments in [39] , since it is not stated there in the precise form that we need. 
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Proof. Theorem 16.11 is essentially Theorem 4.3 from [39] (or rather Theorem 5.1, which deals 
with the periodic boundary conditions instead of no flow - but Theorem 5.1 is a direct corollary 
of Theorem 4.3). We re-scaled time compared to Theorem 4.3 from [39] to make (16.3p hold. 
We also replaced the conclusion of the e-scale mixing (dehned in [39]) with how it is actually 
proved: fl6.4p follows directly from (4.8) and the next estimate in [39] as well as the choice of 6 
immediately below these two estimates. □ 


Here is the key corollary of Theorem 16.11 that we will use in our proof. 

Corollary 6.2. Let /q G C'°°(T^) be a mean zero function. For every e > 0 there exists a Yao- 
Zlatos flow u{x,t) given by Theorem 16'. 11 such that ||Vm||loo ^ 1 for every t and the solution 
f{x,t) of the equation fl6.2p satisfies 

(6-5) ll/(L'?')llil-i ^ C'sll/olU-e 

for some 

(6.6) r ^ C'2e“^'^^| loge|^/^. 

Here C 2,3 ^ 1 are universal constants. 


Proof. To derive this corollary from Theorem 16.11 let us set k = e. We need to address a couple 

and H~^ norm of the solution. 


of issues. The hrst one is the connection between 
Lemma 6.3. Let f G be mean zero. Fix e > 0 and suppose that 

1 


f{x)dx 


< e 


\Qnij\ JQ^ij 

for some n ^ [\ log 2 e|] and i,j = —2"'“^,..., 2”“^ — 1. Then 

(6.7) II/IIh-1 ^C'3||/I|l^6. 


Proof. The proof is by duality. Take any g G H^. Since / is mean zero, without loss of 
generality we can assume that g is also mean zero. Then, denoting gq^ij the average value of g 
over Qnij, we have 


[ fgdx = ^ [ fgdx ^ / f{x){g{x) - gq^.fidx 

i j ^ Qnij i j Qnij 


+ 


y^.gQn 


f{x) dx 




Qr\ 




\\f\\L‘^{Qnij)\\g - gQnij\\L‘i(Qr.ij)^ 




E ||/|U.(0.„)||V9||i.(0..,) + e|l/lll»ll9llL. 




^ C'e||/||L2||c/||^i +e|l/||L° 


m<Ce\\fh^ 


m- 


Here we used Poincare inequality in the last and in the penultimate step, and Cauchy-Schwartz 
inequality in the last step. This proves the lemma. □ 
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Another technical aspect we need to discuss is the smoothness of u. The construction in 
[55] does not explicitly control higher order derivatives of u beyond the Lipschitz condition 
||Vm||loo ^ 1 . However, it is not difficult to see that a properly mollihed velocity held will 
have the same mixing properties up to renormalization by some universal constant. Let 0 be a 
molliher, 0 ^ 0, 0 G C°°(T^), supp(0) C i?i/ 4 ( 0 ), Jj 2 4>ix) dx = 1. Denote (psix) = 6 ~'^(f){x/ 6 ), 
and us{x) = (j )5 * u{x). 

Lemma 6.4. Suppose that an incompressible vector field u{x,t) satisfies ||Vm||loo ^ D for all 
x,t. Assume /o G and denote f{x,t) and fsix,t) solutions of the transport equation 

H 6 . A) with velocity u and mollified velocity us respectively. Fix any T > 0. Then as 6 ^ 0, we 
have \\f{x,t) — fs{x,t)\\L 2 —)■ 0 uniformly in t E [0,T]. 


Proof. The proof of this lemma is standard and elementary. We provide a brief sketch. First 
note that ||VM 5 (-,f)|| 1^00 < ||Vn(-,f)|| ^oo . Consider <Ft(x) and ^t, 5 {x), the trajectory maps 
corresponding to u and us- A straightforward estimate based on Gronwall lemma gives 

Reversing time, we hnd that the same holds for the inverse maps: 

So we get 



\f{x,t) - fs{x,t)\‘^dx 



\foi^f\x))-M^fjix))\^dx ^ 5 ||V/o||i.ce*^. 


□ 


To complete the proof of the corollary, note now that we can choose 6 = (j(/o) small enough 
so that 

\\f5{-,r) - /(•,r)||^_i ^ \\fs{-,T) - /(•,r)||L 2 ^ C'sH/olU^e. 

Then if u{x,t) is Yao-Zlatos vector held yielding fl6.5l) . we can take us as our smooth how and 
get that fl6.5l) holds with the renormalized constant 2 C 3 . □ 


Before stating our main result on Yao-Zlatos hows, we need one more auxiliary result. In 
our scheme, it is convenient to work with the norm of the solution. However Corollary 16.21 
involves the L°° norm, so we need some control over it. We could get it from fl2.6p and Sobolev 
embedding. However the bound in fl2.6p involves norms of higher order derivatives of u and 
would result in weaker estimates for the how intensity needed to suppress blow up. We prefer 
to estimate the L°° norm of the solution directly. 

Let p{x,t) be the solution of fll.21) : 

dtp+{u-V)p-Ap + V ■ p)) =0, p{x,0) = po{x), xeT^. 

where u{x,t) is smooth and incompressible. 

Proposition 6.5. Let 0 ^ p G C°°(T^). Suppose that \\p{-,t) — p \\ l ^ ^ 2B for all t G [0,T] 
and some B 1. Then we also have ||p(-,t) — p\\l°° ^ C4B iaax{B, p^A'j some universal 
constant C 4 and all time t G [0,T]. 


We postpone the proof of this proposition to the appendix. 
We are now ready to state the main theorem of this section. 
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Theorem 6.6. Let po ^ 0 G C'°°(T^), and suppose ||po~pIIl 2 < B for some B > 1. Then there 
exists smooth incompressible flow u{x,t) with ||Vm(- ^ A{B,p) such that the solution 

p{x, t) of the equation (11.21) is globally regular. Here we can choose 

(6.8) A = Cexp (C(l + B + log(l + B + 

for some universal constant C. 

The flow u{x,t) can be represented as 

(6.9) u{x,t) = '^Auflx,t)xiflt), 

j 

where Ij are disjoint time intervals, and Uj are Yao-Zlatos flows given by Corollary\6fM with a 
certain e = e{B,p) > 0 and certain initial data. 

Remark 6.1. We have to deploy different Yao-Zlatos flows in fl6.9l) dne to the fact that these 
flows are designed to mix a specific initial data, and one can envision nonlinear dynamics 
attempting to break the barrier in different ways. 

Remark 6.2. Similarly to Theorem 15.41 for the RE flows, one can nse combinations of Yao- 
Zlatos flows to achieve stronger results, such as convergence of the solution to the mean, and 
at increasingly fast rate if the flow amplitude is allowed to grow. We will not pursue a walk 
through these results here, instead leaving it to the interested reader. 


Proof. The scheme of the proof is similar to the RE flows case, but Corollary 16.21 replaces 
Lemma [5.31 with some necessary adjustments. 

We start with n = 0 on some initial time interval. If ||p(-,t) — p||l2 ^ B for all t, global 
regularity follows. Suppose this is not the case, and let to be the first time when ||p(-, to)—p\\L 2 = 
B. Similarly to the RE case, we also know that ||p(-, ^o)||pi ^ Bi. In addition. Proposition 16.51 
ensures that ||p(-,to) — p|| ^ C 4 B max(i?, p^'^^). 

Fix e given by 

(6.10) e =- -. 

SCs^ACoB^ + 2p + 1(1 + C 4 max(R, pV2)) 

Take a flow u guaranteed by Corollary 16.21 corresponding to this value of e and the initial density 
p{x,to), and let r be the time in (16.6j) . Denote ri^{x,t) the solution of the equation 

( 6 . 11 ) dtp^ + A{u-V)r]^ = 0, p^{x,0) = po = p{x,to). 


Then by Corollary 16.21 and Proposition 16.51 we have 


( 6 . 12 ) 


- p\\h-i ^ CsWpo - 






B 

8-\/ ACqB^ 2p -1- 1 


We will set u{x, t) = 0 in (I6.1ip for t ^ t/A. 

Now we are going to turn on the same flow u{x, t) in the equation for p^ at time tQ, for the 
duration r/A. Let us denote p^{x,to + t) the solution of the equation 


a^p^-F R(m ■ V)p^ - Ap^-h V ■ (p^V(-A) ^(p^-p)) = 0, p^(x,0) = p(x,to)- 
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The first condition that we are going to impose on A is that 

2r 


(6.13) 


A 


^ Cl min (l, p ^,B . 


Given fl6.6D and fid.iop . it is easy to check that (I6.13P holds if fl6.8p is satished. Recall that 
given fl6.13l) . Proposition 13.11 ensnres 

(6.14) jjp^(x,to+ i) - Ph^ ^ 2B 

for t G [0, 2r/R] and the solntion stays smooth on this time interval. Next, by the approximation 
Lemma [4.21 we have 

(6.15) 

^l|p"^(-Go+ ^) -h^(-G)lli 2 ^ -||p^(-Go+ ^)||^1 +4||p(-,to)||^iexp \\\/u\\L^ds 

+ C'||p"^(-, to + t) — p||^2 (||p"^(-, to + t) — PIII 2 + p^) 

^ ABl exp{2CAt) + ACB^ {{2Bf + p^) 

for every t G [0,2r/R]. Here we nsed fl6.14p in the last step. Let us now impose the second 
condition on A which says that it should be large enough so that 

(6.16) ^exp(4Cr) + ^4013^ ({2Bf + P) < 

Note that in particular (I6.16p implies that 

(6.17) 0.8R ^ ||p^(-, h + t)- pWl^ ^ 1.2B 


for t G [0,2r/H]. Also, fl6.16p . fl6.12p and fl6.15p can be used to estimate that for every t G 
[t/A,2t/A] we have 

\\p^{-Ao + t) -p\\h-i ^ llh^(-G) -pllit-i + l|p^(-Go + t) -h^(-G)IU2 ^ 

B B _ B 

8^y4:CoB^ + 2p+l 8V4GoH2 + 2p + l 4^4GoH2 + 2p + 1 ‘ 

Therefore, using fl6.17p . we obtain 

(6.18) ||/(., f„ + f) - (ill*. ^ 3^ 2BV4C„B^ + 2p+l 

IIP f-,to + t) — Pllij-1 

for all t G [t/A,2t/A]. The proof of global regularity is completed similarly to the proof of 
Theorem 15.21 using fl3.4p . (16.17p and fl6.18p . 

Finally, the sufficiency of the condition (16.8p follows from straightforward analysis of the 
bounds (I6.10p . (16.6p . and (I6.16p . The approximation lemma constraint (I6.16P is what truly 
determines the exponential form of (16.8p . □ 


Remark 6.3. One can check that by using Theorem 5.4 from [39] and similar arguments, we 
can change the periodic setting to the hnite domain with no-slip or no-flow boundary condition 
for u{x,t), and obtain analogous results but with a somewhat weaker estimate for the flow 
intensity. We leave the details to the interested reader. 
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7. Discussion and generalization 


The scheme developed in the previous sections of this paper should be flexible enough to be 
applied in different situations. Here we briefly and informally discuss the main features that 
appear to be necessary to apply our analysis. On the most general informal level, one can say 
that the idea of the scheme is that the fluid flow, if sufficiently intense and with strong mixing 
properties, can make a supercritical equation into subcritical one for a given initial data. 

It appears that for the scheme to be applicable to a nonlinearity N{p) we need that for the 
solutions of the equation 

dtp + {u ■ V)p - Ap + N{p) = 0, 

either the mean or some norm of p does not grow or at least obeys global hnite (even if growing) 
bound in time. Without such assumption, it is difficult to rule our hnite time blow up of the 
mean value of the solution, which large norm has no way to arrest. The second condition 
that is needed concerns the bound on the nonlinear term in the spirit of 


(7.1) 


j N{p)pdx 


<cf{\\p\\L^)\\pnt, 


with a < 2. This would allow control of the norm growth by diffusion when norm is 
large. The third condition is that some analog of the approximation lemma holds. This seems 
to require bounds similar to fl7.1l) . 

Of course, the scheme can also be adapted to the cases where diffusion is given by some 
dissipative operator other than Laplacian, for example a sufficiently strong fractional Laplacian, 
in which case the norm needs to be replaced by some other norm natural in the given context. 
It is also likely that diffusion term does not have to be linear, even though this may require 
subtler analysis. 

As far as other possible classes of flows that may have the blow up arresting property, the 
main clearly sufficient requirement for our scheme to be applicable appears to be as follows. 
First, the flows should be sufficiently regular and in particular satisfy Lipschitz bound in spatial 
variables. Secondly, for every e > 0 it should be possible to hud a flow from the given class, 
with uniform in time Lipschitz bound, such that for every /q G C°° the solution f{x,t) of the 
transport equation 

dtf + (we ■ V)/ = 0, f{x, 0) = fo{x) 

satishes 


(7.2) ll/(•,A)||^_7 ^eC(||/o|U») 

for some < oo. Observe that even though we did not frame the discussion of the mixing 
effect of the RE flows in terms of the norm. Lemma 15.31 clearly implies that fl7.2p holds 
for the RE flows. There are other classes of flows that look likely to satisfy these properties, 
such as optimal mixing flows discussed in [25]. Generally, decay of the H~^ norm is one of the 
general measures of the mixing ability of the flow, hence we have a link between efficient mixing 
and suppression of blow up, which is quite natural. We refer to m, m, m, imi for further 
discussion of H~^ norm as a measure of mixing and some bounds on mixing rates for natural 
classes of flows. It also looks possible that some flows that do not in general lead to norm 
decay without diffusion can still be effective suppressors of blow up if diffusion is taken into 
account. A natural and common class to be investigated here are some families of stationary 
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cellular flows. Furthermore, similarly to [5] , our construction can be applied more generally to 
the case where the transport part of the equation is replaced by some other unitary evolution 
for which an analog of (I7.2jl holds. We plan to address some of these generalizations in future 
work. 


8. Appendix I: Finite time blow up 


In this section, our main result is a construction of examples where solutions to the Keller- 
Segel equation set on without advection fll.ip blow up in hnite time. As we mentioned in 
the introduction, similar results are well known in slightly different settings. The argument 
below is included for the sake of completeness. It is closely related to the construction of [33], 
but is simpler. The argument is essentially local and can be adapted to other situations as well. 
Although we will focus on the d = 2 case, some auxiliary results that remain valid in every 
dimension will be presented in more generality. 

Theorem 8.1. There exist pQ G C°°(T'^), po ^ 0, such that the corresponding solution p{x,t) 
of the equation M.l\) set on blows up in finite time. 


Without loss of generality, we will assume that the spatial period of initial data and solution 
is equal to one, so = [—1/2,1.2]^. Let us hrst state a lemma that will allow us to conveniently 
estimate the chemotactic term in the equation. 


Lemma 8.2. Assume = [—1/2,1/2]'^, d^ 2. For every f{x) G we have 


( 8 , 1 ) 


V(-A)-‘(/(U-/) =-lim 

Cd 1-*"+, 


(x - p) 

\x — y\'^ 


{f{y) - /)e 


Here on the right hand side f{y) is extended periodically to all / denotes the mean value 
of f, and Cd is the area of the unit sphere in d dimensions. 


The expression fIS.ip is of course valid for a broader class of /, but the stated result is sufficient 
for our purpose. 


Proof. Without loss of generality, let us assume that / is mean zero. By dehnition and properties 
of Fourier transformation, we have 


V(-A)-V(x) 



^2'Kikx 




ik 

2'x\k\^ 


/>)• 


To link this expression with fIS.ip . observe hrst that for a smooth /, a straightforward compu¬ 
tation shows that 


( 8 . 2 ) 


- E » 




= — lim 

7^0+ 


^27vipx 


ip 

27r|pP 


e-^^^Py-^\y\"fiy)dydp, 


where the function f{y) is extended periodically to the whole plane. Indeed, all we need to 
do is plug in the Fourier expansion f{y) = integrate in y, and observe that 

(Tr/y)'^/^ exp(—7r^|/c — p\/'j) is an approximation of identity. 
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On the other hand, recall that the inverse Laplacian (—A) of a sufficiently regular and 
rapidly decaying function g is given by 


(8,3) 


^2Ttipx 


e-^^^^yg{y)dydp = 


lRd^og\x-y\g{y)dy d = 2; 
j^J^d\x-y\‘^~'^giy)dy d ^ 3. 


jRd (27r|p|)2 

The expression on the right hand side of (18.21) . with help of (I8.3p . can be written as 


Right hand side of (18.211 = lim 

7->0+ 

Integrating by parts, we obtain (18.ip . 


^ /]R 2 log 1^ - dl V (^f{y)e dy 


d = 2; 
d ^ 3. 


□ 


Suppose that the initial data po is concentrated in a small ball Ba of radius a centered at 
the origin, so that Jg podx = Jj 2 Podx = M. Suppose that ^ > fe > 2a, M > 1, and let 0 
be a cut-off function on scale b. Namely, assume that 0 G 0°°(T^), 1 ^ 0(x) ^ 0, 0 = 1 on 
Bb, and 0 = 0 on The function 0 can be chosen so that for any multi-index a G Z^, 
|R"0| ^ The parameters a, M and b will be chosen below. The local existence of 

smooth solution p{x, t) can be proved by standard method, see e.g. [22] for a closely related 
argument in the setting. It is straightforward to check using parabolic comparison principles 
that if Po ^ 0, then p{x,t) ^ 0 for all f ^ 0. Also, we have Jj .2 p(x,t) dx = M, at least while 
p(x, t) remain s smooth. 

The hrst quantity we would like to consider is J ^2 p{x,t)(j){2x) dx. We need an estimate 
showing that the mass cannot leave Bb too quickly, we note that the constants Ck employed 
later in this section are not related to the constants Ck in the previous sections. 


Lemma 8.3. Suppose that a, b, 0, M and po are as described above. Assume that the local 
solution p{x,t) exists and remains regular in the time interval [0,r]. Then we have 


(8.4) 

for every t G [0, r]. 

Naturally, the bound 
Proof. We have 


dt / p{x)(p{2x) dx = 


'T2 


p(x, t)0(2x) dx ^ M — CiM^b ^t 


is only interesting if t is sufficiently small. 


'T2 


Ap{x)(f){2x) dx 


Ij2 


(j){2x)V ■ {p{x)V{—A) {p{x)—p))dx. 


First, using the periodic boundary conditions and integrating by parts, we hnd that 


(8.5) 


Ij2 


Ap{x)(j){2x) dx 


= 4 


It2 


p{x)A(j){2x) dx 


< CMb 


-2 


Next, let -0 G be a cutoff function, if{x) = 1 if |(r| 0 1/2, 0(x) = 0 if |x| ^ 1, 

0 ^ 0((r) 0 1, |V0((r)| ^ C. Using Lemma I5T21 we have 



0(2a:)V ■ (p((r)V(—A) ^{p{x)—p))dx 
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_ 1 
TT 

^ c 
+ c 


/ '\/(p{2x)p{x,t) lim 

Ij2 7-^0+ ^ 


jx-y) 

\X - |/|2 


(piVyi) - p)e dydx 


/T2 


V(j){2x)p{x) 


/T2 


jx-y) 

lx — I/P 


(piv) - p)i^{y) dydx 


V(j){2x)p{x) lim f j^3^(p(l/) - p)e ^'^''(1 - i’iy)) dydx 


/ l’2 7—>-0 + 

= C{I) + C{II). 

We passed to the limit 7 —)■ 0+ in the hrst integral since the integral of the limit converges 
absolntely. Using symmetrization, we can estimate 


iI)<P 
+ 



(V0)(2x)p(x,t) 
p{x)p{y) 


{x - y) 


'il^{y) dydx 


x-y)- {V(j){2x)i/j{y) - V(p{2y)t/j{x)) 


dxdy 


7 CM%-^ + 
7 CM%-^. 


|x - I/P 

P(a;7)p(2/7)(||VVIU- + ||V0||loc||V'0||loo) dxd|/ 




We use the fact that supp0 C suppi/^ C Bi in the second step. 

Next let us estimate {II)- Note that in this case the kernel is not singular, since supp0 C i?2fe 
while supp(l —'ijj)<Z B^. However, there is an issue of convergence of y integral over the inhnite 
region. Suppose that x G i?26- Dehne mean zero function g via p{y) — p = Ag. In fact, we 
have g{k) = —/3(/c)(27r|/cp)“^ for k ^ 0. By working on the Fourier side, it is easy to show that 
||fi'IUi(T2) 7 ||5 'I|l 2 (t 2 ) 7 C'||p|7i(t2). Now we can estimate 

[ W— ^{Piy^t) - p)e-^^y^\l-^|J{y))dxdy 

1-^0+V \x-yr 

= lim [ ^ — ^Ag{y,t)e-'^^y^\l-'ijj{y))dxdy 
^^0+iR2 |X-|/P 

+ 2V (jTipp) ^ - ’/'(»))) + |P_ ^2 ^ - </)(»))) j dxdy. 

Here g{y,t) is extended periodically to the whole R^. Note that in the hrst summand in the 
last integral in (18.dh we can pass to the limit as 7 —)■ 0 since A ^ \xly \2 j decays sufficiently fast. 
For every x G Bb, we obtain an integral which is bounded by U|||/|7i |n|>o ^ CM. 

It is straightforward to estimate that the last two summands in fl8.6p are bounded by 


C [ \g{y,t)\{j\y\ ^ + 7^||/|)e -'^(l/)) d|/7 
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C|kllLi(T2) ^+ 7Vl)e ^ 0. 

n€Z^,|n|>0 

Combining these estimates, we see that 

(JJ) ^ CM [ {\/(p){2x)p{x,t) dx ^ CM%-\ 

Therefore, for all times where smooth solution is still dehned, and under our assumptions on 
values of parameters, we have 


dt 



p{x, t)(l){2x) dx 


^ CM%-\ 


This implies (18.dh and hnishes the proof of the lemma. 


□ 


Let us now consider the second moment J ^2 dx. Closely related quantities are 

well-known tools to establish hnite time blow up in Keller-Segel equation; see e.g. [35], |33] . 
We have the following lemma. 

Lemma 8.4. Suppose 1/4 ^ 5 > 0 and (j) is a cutoff function on scale b as described above. 
Let pq e and assume that the unique local smooth solution p{x,t) to U.l\) set on 

is defined on [0,T]. Then for every t G [0,T] we have 

{8.7) dt J \x\‘^p{x,t)4>{x)dx ^p{x)4>{x) dx'^ + C 2 M\\p\\Li(j 2 \Bt) + C^bM'^ + C^M. 

Proof. In the estimate below, we will use the formula (18.Ih with 7 set to be zero. All the 
estimates can be done completely rigorously similar to the proof of Lemma [H31 we will proceed 
with the formal computation to reduce repetitive technicalities. 

We have 


dt / \x\‘^p{x,t)(p{x) dx = / \x\‘^Ap{x)(p{x)dx 


/T 2 


/T 2 


+ / \xf(j)V ■{pV{-A)-fip-p))dx 

JT 2 

= 4 / <f)pdx+ / \x\^A(f)pdx + A / {x-Vffipdx 


If 2 


If2 


- / fi{x)p{x) 
'X Jf2 

7^ [ \x\‘^Pix) 


JT 2 

— — ^(piv) - p) dydx 


271 


/T 2 


|x - y\ 

V(p{x) -jx-y) 
|x -1/|2 


{p{y) - p) dxdy. 


= (i) + (a) + (Hi). 
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Here (z) denotes the first three terms. By our choice of (p, {i) does not exceed C 4 M for some 
constant C 4 . Next, let us write 

(ii)=--f (p(x)p{x,t) [ ^- p)^p{y)dydx 

/Q Q^ ^ ~ y\ 

lO-OJ . n 

-/ Hx)p{x,t) ^ — ^{p{y,t) - p){l-ilj{y))dydx, 

'X Jr2 |x — I/p 

where "0 is a cutoff function as in Lemma 18.31 The absolute value of the integral 
/ (l>{x)p{x,t) f ^ — ^{p{y,t) - p){l-'il){y))dydx 

Jt^ Jir 2 \x-yr 

can be controlled similarly to the estimates applied in bounding the term {II) in the proof of 
Lemma 15751 leading to an upper bound by CM%. Next, we can estimate 


P [ ^{x)p{x,t) [ ^ — ^^p{y)dxdy 
JT2 JR2 \x — I/p 


as well. Split the remaining part of the hrst integral in 


^ CM^b 

into two parts: 


(l>{x)p{x,t) [ ^ — ^p{y,t)(p{y)dxdy 
7rJT2 JRS \x-yr 

[ 4>ix)p{x,t) [ ^ — ^p{y,t){l-(j){y))^l){y)dxdy 
71 Jt2 jr2 f - i/r 


Using symmetrization, we obtain 
1 


/ 

(t){x)p{x) 

Jt2 


1 

/ / 


x{x — y) ,,,,,, , 

- wpiy)^iy) dydx 


|x - I/p 


1 


'T2 


p{x)(j){x) dx 


On the other hand. 


/T 2 


Hx)p{x,t) I f)(l - (t>{y))'^{y) dxdy 



'T 2 

p{x,t)p{y,t) 

\x — y\^ 


x-y)- [x(j){x){l - (p{y))i>{y) - y4>{y){l - (p{x))i>{x)] dxdy. 


Let us dehne F{x,y) = x(f){x){l — 0(|/))|p(|/) — y(f){y){l — (f){x))'ip{x). Observe that F{x,y) = 0 
on Bb X Bb, F{x,x) = 0 and \VF{x,y)\ ^ C for all x,y. This means 

\F{x,y)\ = \F{x,y) - F{x,x)\ ^ \\VF\\l^\x - y\xBixBi\B,xB,{x,y) 

^ C\x - y\XBixBi\Bi,xBf,{x,y), 

where Xs{x,y) denotes the characteristic function of a set S' C x R^. Therefore, 

x{x — y) 


Ij2 


(t>{x)p{x,t) 


\x-y\ 


-p{y)0--(t>{y))'fp{y)dxdy 
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p{x,t)p{y,t)dxdy 
J J 

To summarize, {ii) can be bounded above by 

p{x,t)(j){x) dx^ + CM\\p{-,t)\\Li(j 2 \B,) + CbM"^. 


1 


'T2 


Finally, let us estimate (Hi). Similarly to the previous part, we have 


112 


\x\‘^p{x,t) [ {p{y, t) - p)(l - ^P{y)) dxdy ^ CbM^. 


Also, 


P / \x\‘^pix,t) 
Jt 2 


\x-y\ 
f \/(p{x) -{x-y) 


(1 — 'ip{y)) dxdy ^ CbM'^ 


\x — y\^ 

as well. The remaining part of {Hi) we can estimate by using symmetrization; 

V(j){x) ■ {x — y) 


\xpp{x, t) 



p{x,t)p{y,t) 


12 -p{y^t)ij{y)dxdy 

X - y\^ 

x-y)- (V0(x)V'(|/)|xp - V(t){y)'il^{x)\y\^) 


\x - y\^ 


dxdy. 


Observe that 

\\x\‘^V(j){x)'ijj{y) - \y\‘^V(l){y)'ip{x)\ ^ CxB2,xB2b\B,xBjx - y\. 

Therefore 

{in) ^ CbM"^ + CM\\p\\ii(j 2 \B,,)- 

Combining the estimate of {i), {ii) and {Hi) yields (18.7p . proving the lemma. □ 

We are now ready to complete the proof of Theorem 18.11 

Proof of Theorem \8.1[ Let us recall that we assume 1/4^6^ 2a, and the initial data po is 
supported inside Ba- Assume that the unique solution p{x, t) of fll.ip set on remains smooth 
for all t. Then by Lemma 18.31 and conservation of mass, for all t ^ 0 we have 


||p(-,t)||Li(T 2 \ij,) ^ M - / p{x,t)(j){2x)dx ^ CiM^b H. 


'T2 


Also, 


/T2 


p{x,t)4>{x)dx ^ / p{x,t)4>{‘2x) dx ^ M — CiM^b 


/T2 


Therefore, by Lemma 


we have that 


dt I \x\‘‘^p{x,t)(j){x) dx ^ --^{M — CiM'^b 


/T2 


27r 


+C3M% + CiM 


for all 0 ^ ^ 


b2 


CiM 

1. Choose b so that C^b ^ 0.001. 


We will now make the choice of all our parameters. 
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2. Choose M so that M ^ IOOOC4. 

3. Choose a so that the following three inequalities hold: 

a ^ 6 / 2 , a ^; - , and a ^-^-=. 

4. Choose the time r = 

With such choice of parameters, it is straightforward to check that 

r ]\^2 

(8.9) dt \x\^p{x,t)(l){x)dx ^—— 

50 

for every t E [0, r]. But by assumption, supppo C Ba, and so 

( 8 . 10 ) / \x\^P q{x)(1){x) dx ^ a^M. 

Together, fl8.9p . fl8.10p and our choice of r imply that \x\‘^p{x,T)(j){x) dx must be negative. 
This is a contradiction with the assumption that p{x,t) stays smooth throughout [0,r]. □ 

In fact, it is not hard to verify that hnite time blow up persists if we add an advection term 
to the Keller-Segel equation, but change the order of the selection of u and po compared with 
the results on suppression of chemotactic explosion. Namely, the following theorem holds. 

Theorem 8.5. Consider the Keller-Segel equation (II.2p set on T^. Suppose that u{x,t) £ 
C°°(T^ X [0, 00 )) is incompressible. Then there exists po G C°°(T^), po ^ 0, such that the 
corresponding solution p(x, t) of (II.2p blows up in finite time. 

The proof of Theorem 18.51 closely follows that of Theorem 18.11 The advection term can be 
easily estimated and the proof requires only a few adjustments of constants. We omit the 
details. 


9. Appendix II: Some inequalities 

Here we hrst prove Proposition 16.51 and then sketch the proofs of the multiplicative inequal¬ 
ities ( 12 . 2 p and ( 12 .5p . 

Let p{x, t) be a solution of (11.21) 

dtp + (u ■ V)p - Ap + V ■ (pV(-A)“^(p - p)) = 0, p{x, 0) = po{x), 
with some smooth incompressible vector held u{x,t). 

Proposition 9.1. Let po G C°°(T^). Suppose that ||p(-,t) — p\\l^ ^ 2B for all t G [0,T] and 
some B ^ 1. Then we also have ||p(-, t) —p||l°° ^ C 4 -B max(i?, p^/^) for some universal constant 
C 4 and all t G [0, T], 

Proof. Observe hrst that by a direct computation, for every integer p ^ 1 we have 

dt [ (p-pf^dx =-{A--) [ |V((p-p)P)|^da; 

+ {2p — 1) / (p — pY^~^^ dx -I- 2pp f (p — pfi^ dx. 

JT2 Jt2 


(9.1) 
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To obtain fl9.ip . we need to mnltiply fll.2p by 2p{p — integrate, and then simplify the 

obtained terms using integration by parts and the fact that u is divergence free. 

Let us estimate ||p — p\\l 2 ^ inductively. By assumption, we have ||p(-,t) — p\\l^ ^ 2B for all 
t G [0,T]. Assume that for some n ^ 1, we have 

IIp-pIIl 2" ^ for alH G [0,T]. 

Let us derive an estimate for an upper bound T„+i on \\p — p\\ 12 ^+^ on [0,T]. For that purpose, 
let us set p = 2” in fl9.ip and let us dehne /(x, t) = [p — py = (p — Then fl9.1l) implies 

(9.2) dt f \f\^dx^-2[ \Vf\^dx + 2^+^ [ |/|2+2-"dx + 2"+ip / I/p. 

JT2 JT2 JT2 JT2 

Also, in terms of /, our induction assumption is that J ^2 \ f\dx ^ 

We will now need the following Gagliardo-Nirenberg inequality. 


Lemma 9.2. Suppose v G d ^ 2, and the set where v vanishes is nonempty. Assume 

that q,r > 0, 00 > q > r, and ^ | > 0. Then 


(9.3) 


v\\lo ^ G(d,g)||Vn||/2 


l—a 
L’’ ) 


a = 


1 _ 1 

r q 


1 

d 



The constant C{d,q) for a fixed d is hounded uniformly when q varies in any compact set in 
(0,cx)). 


Proof. This inequality is well known in the case v G C'^(M'^), see e.g. |ni]. A simple proof is 
contained in ra- Going through the proof in |22], it is not difficult to verify that the result 
still holds in the periodic case under the assumption that v vanishes somewhere in T'^ (which 
rules out increasing mean value without increasing variance). One can similarly trace the claim 
regarding the constant C{d,q). We refer to [22] for details. □ 


Applying Lemma [9.21 with d = 2, r = 2, and g = 2 + 2 ” yields 

(9.4) ll/fpiA $ C||V/||ir||/||i. ^ ^11 V/lli. + CWfWif^, 

where we used Young’s inequality in the last step. Moreover, we also have 

(9.5) ll/IU-$C||V/|ll^||/|liy 
Applying fl9.5p and fl9.4p to fl9.2p . we obtain 

(9.6) S,||/||i, s: -C-,||/||L||/|l2? + C,2”+‘||/||ir^ + 2”+y||/||i=, 

where Ci ^2 are some hxed universal constants (not connected to Ci and C 2 used earlier in the 
paper). Glearly, given the upper bound on ||/||li, the right hand side of (19.6p turns negative 
if ||/||l2 becomes sufficiently large. Thus ||/||l2 cannot cross this threshold. Assuming without 
loss of generality that ^ 1 for all n, a direct computation shows that if ||p —p||£2"+i reaches 
the value T„+i which satishes the following recursive equality, then the right hand side of fl9.6p 
is negative: 


logT„+i = max(r„,0„) 
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where 

2*1+1 _ ^ 

(9-7) ^logT„ + ^((n+ l)log2 + logC') 


(9.8) Qn = logT„ + ^^((n + l)log2 + logC + max(logp,0)). 

Here C ^ 1 is some universal constant. Denote qj = and observe that due to telescoping, 


n 




2*1+1 _ 
2*1 


n—>-oo 


> 2 . 


An elementary inductive computation shows that if H then the hrst recursive relation 

fl9.7p determines the size of T„+i, yielding the estimate T„+i ^ CB^. If B < then the 
second relation fl9.8p dominates, yielding the estimate T„+i ^ CBp^^'^. Since 

IIp-pIIl“ = lim || p - p || l 2 * i , 

n—>-oo 

we obtain that 

Up - p\\l°° ^ CB max(i?, p), 

proving the proposition. □ 


Proposition 9.3. Suppose that f E C'°°(T“) and is mean zero. Then 


1—a 






a = 


^-j + l 

n 


where D stands for any partial derivative, 2 ^ p ^ oo, and we assume n > m + d/2. 


Note that the last assumption is not necessary. However it makes the proof simpler, and 
this is the only case we need in this paper. Indeed, in the estimates of Section 3 we have 
s + 1 > I + d/2 unless I = s (recall d ^ 3). But when I = s, we are only estimating ||D^p||2,2, 
which is straightforward. 


Proof. Consider p = 2. Then 

-t m m 

(9.9) p™/iu^^ii/iil;"ii/ii^. 

by Holder inequality on Fourier side. 

Next consider p = oo. Then 

53 |A^ri/+)i+c 53 iiri/+)i = (/) + (//). 

0<|fc|<A \k\^A 


Now 


(/) ^ C'A”^+5 


53 mw 

0<|fc|<A 


1/2 
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by Cauchy-Schwartz. On the other hand, 

(//) < c |4ri/(4)|A s;c||/||«„a‘’"-”>A, 


\k\-^K 


JfcIM 


provided that n > m + |. Choose A so that 


L2 


A "'+2 = 


H' 




Snch choice leads to the bonnd 

n — m-\-d/2 m+d/2 

(9.10) WD-fh^ ^ C^II/IIZ^^II/li;:^- 

The general case 2 < p < oo follows immediately from flQ.Op and flO.lOp . 


□ 


Proposition 9.4. Suppose that f E and m > 0. Then 

Proof. The proof of this proposition can be done similarly to the previous one. One needs to 
use that ||/||l°° ^ ll/IUi- We leave details to the interested reader. □ 
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